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Abstract 

Two variable flavor number schemes are used to describe bottom quark pro- 
duction in deep inelastic electron-proton scattering. In these schemes the 
coefficient functions are derived from mass factorization of the heavy quark 
coefficient functions presented in a fixed flavor number scheme. Also one has 
to construct a parton density set with five light flavors (u,d,s,c,b) out of a set 
which only contains four light flavors (u,d,s,c). In order a 2 the two sets are 
discontinuous at /i = nib which follows from mass factorization of the heavy 
quark coefficient functions when it is carried out in the MS-scheme. Both vari- 
able flavor number schemes give almost identical predictions for the bottom 
structure functions i*2 & an d Fl Also they both agree well with the corre- 
sponding results based on fixed order four-flavor perturbation theory over a 
wide range in x and Q 2 . 

PACS ll.lOJj, 12.38Bx, 13.60Hb, 13.87Ce. 
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I. INTRODUCTION 



The HI and ZEUS experiments at HERA now have enough integrated luminosity to 
study bottom quark deep inelastic electroproduction. Therefore it is an appropriate time 
to present up-to-date predictions for the bottom quark components of the deep inelastic 
structure functions Fi(x, Q 2 ) where i = 2, L. We assume that the bottom quark is produced 
in an extrinsic fashion and that the neutral current reaction dominates over the charged 
current one. This means that in fixed order perturbative QCD the heavy quark structure 
functions F ijb (x, Q 2 , m|), i = 2, L are given by the virtual-photon gluon-fusion processes and 
their higher order corrections with only light partons in the initial state. Notice that in the 
case of bottom quark production the light partons are represented by the gluon and the four 
light flavors u, d, s, c together with their anti-particles. 

In the literature one has adopted two different treatments of extrinsic bottom quark 
production, which are known as the massive and massless descriptions. The former treats 
the bottom quark as a heavy quark (with mass m&) and the partonic cross sections (or 
heavy quark coefficient functions) are described by fixed order perturbation theory (FOPT) 
as mentioned above. Notice that due to the work in ffl] the perturbation series is now known 
up to second order. The latter treatment, which has been rather popular among groups 
which fit parton densities to experimental data, treats the bottom quark as a massless 
quark so that it can be represented by a scale dependent parton density fb(x, /x 2 ). Although 
at first sight these approaches are completely different they are actually intimately related. 
It was shown in || that the large logarithms of the type ln(Q 2 /m 2 ), which appear in FOPT 
when Q 2 ^> m 2 , can be resummed in all orders. The upshot of this procedure is that 
the bottom components of the deep inelastic structure functions F ib (x,Q 2 ,ml), where i = 
2, L, which in the FOPT approach are written as convolutions of heavy quark coefficient 
functions with four-flavor light-mass (u,d,s,c) parton densities, become, after resummation, 
convolutions of light-mass parton coefficient functions with five-flavor light-mass parton 
densities which also include a bottom quark density. This procedure leads to the so-called 
zero mass variable flavor number scheme (ZM-VFNS) for Fi t f,(x, Q 2 ) where the mass of the 
bottom quark is absorbed into the new five flavor densities. To implement this scheme one 
has to be careful to use quantities which are collinearly finite in the limit m b — > 0. From 
the above considerations it is clear that the FOPT approach is better when the bottom 
quark pair is produced near threshold (where Q 2 (x~ 1 — 1) > 4m 2 ) because terms in m& 
are important in this kinematic region. On the other hand far above threshold, where also 
Q 2 ^> ml, the large logarithms mentioned above dominate the structure functions so that 
the ZM-VFNS approach should be more appropriate. Both approaches are characterized by 
the number of active flavors involved in the description of the parton densities which are 
given by four and five respectively. Each scheme has different densities but the momentum 
sum rule either gets contributions from four-flavor densities or five-flavor densities and is 
always satisfied. 

As most of the experimental data will be in the kinematical regime which is between 
the threshold and the large Q 2 region, a third approach, called the variable flavour number 
scheme (VFNS), has been introduced to describe the heavy quark components Fi } b(x, Q 2 ) of 
the deep inelastic structure functions. Actually there are several such schemes. They include 
the Aivasis, Collins, Olness, Tung (ACOT) ||, scheme, the Buza, Matiounine, Smith, van 
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Neerven (BMSN) |2j, |4[] scheme, the Thorne, Roberts (TR) scheme and the Chuvakin, 
Smith, van Neerven (CSN) || scheme. A discussion of them is given in the last reference. 
The difference between the schemes can be attributed to two ingredients entering in their 
construction. The first one is the mass factorization procedure carried out before the large 
logarithms can be resummed. The second one is the matching condition imposed on the 
heavy quark density, which has to vanish in the threshold region of the production process. 
Another aspect of these approaches is that one needs two sets of parton densities. For 
bottom quark production one set only contains densities in a four-flavor number scheme 
whereas the second one, which also includes a bottom quark density, is parametrized in a 
five-flavor number scheme. Both parameterizations have to satisfy the matching relations 
quoted in f2j. Up to next-to-leading order (NLO) they are continuous at the scale /i = mi, 
whereas in next-to-next-to leading order (NNLO) the parton densities become discontinuous 
while going from a four to a five flavor scheme. Starting from a three-flavor set of parton 
densities given in |7J we have recently constructed in a four-flavor set of densities which 
satisfied the matching relations in at the scale // = m c . Then we evolved these densities 
with LO or NLO splitting functions up to the scale /J, — m& and constructed a five-flavor 
set which also satisfied the matching relations in [0. This set was further evolved with 
LO and NLO splitting functions up to high scales. Notice that since the NNLO splitting 
functions are unknown the only difference between the NLO and NNLO parton densities can 
be attributed to the boundary conditions at fi = m c and /i = m b where the latter densities 
become discontinuous contrary to the LO and NLO ones. We can now use these densities to 
discuss VFNS for bottom quark deep inelastic electroproduction, in particular in the CSN 
and BMSN schemes. The previous discussions in were focussed on applications to charm 
quark electroproduction. 

Since any description for bottom quarks follows closely that for charm quarks we refer 
the interested reader to |J for most of the details and simply specialize to bottom quark 
electroproduction in Sec. II. We work to second order in the running coupling constant a s (/x 2 ). 
Numerical results are shown for the structure functions i*2,& and Flj> i n the CSN and BMSN 
schemes. 
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II. BOTTOM QUARK STRUCTURE FUNCTIONS 



In this Section we consider bottom quark deep inelastic electroproduction in two variable 
flavor number schemes, namely the BMSN scheme as proposed in 0, Q and in the CSN 
scheme as proposed in ||. For that purpose we have constructed in || a five-flavor parton 
density set from a four-flavor parton density set. Using our densities we will study the differ- 
ences between the bottom components of the deep inelastic structure functions F^ SN (rij + l) 
and 6 (n/ + 1), where the number of light flavors is nf — 4. 

To keep the discussion short we simply refer the reader to Sec. Ill of for a discussion of 
the MS parton densities, the exact solution for the running coupling constant and the scale 
choice. All references to three-flavour (four-flavor) densities should be replaced by four-flavor 
(five-flavor) densities respectively. All our calculations of next-to-leading (NLO) and next- 
to-next-leading order (NNLO) quantities use A|f 5i6 = 299.4,246,167.7,67.8 MeV, which 
yields a s (5,M|) = 0.114. The structure functions are defined in Eqs.(3.9)-(3.17) of ||, 
where now rij = 4 and m c is replaced by mj = 4.5 GeV/c 2 . To make this paper reasonably 
self-contained we now reproduce the final formulae we use for the structure functions. For 
i = 2, L the CSN scheme uses 

fSf(n / + l,A,Q 2 ,m 2 )=4 
+ a,(n, + 1, {/«$(„, + 1, „*) ® 

+ /r LO (.v + i,^)®c^'"(g,^)j 

W,{n, + 1, {/&,(„, + 1, (C< 2 > (n, + 1, ^ ^) 

+/ 5 s ' LO ("/ + ^ 2 )«C N ' s ' (2, (^,f)[ 

+a*(n / + 1, //) £ el /J&(n, + 1, f) ® C NS ' (2) ^) • 
fe=i m 

where we choose the heavy quark Q to be the bottom quark, so that the number of light 
flavors is rif = 4. The charge of the bottom quark is cq = —1/3 and its mass is m = m&. 
The coefficient function Lf^ RD depends on a parameter A which refers to the invariant 
mass of the Q (Q-pair. For bottom quark production we choose A = 100 (GeV/c) 2 . The 
coefficient functions labelled by C CSN depend on the heavy quark mass but are finite in 
the limit m —>■ 0. They are defined in Eqs.(2.8)-(2.20) in |J. To simplify the notation we 



/^T(«/ + ^ 2 )cS N ' NS ' (0) (^: 
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will refer to the above structure functions as Fp^ N (nf = 5), indicating that they depend on 
five-flavor parton densities. 

The same parameters eg, m b , A etc., also show up in the expressions for the bottom 
quark structure functions in the BMSN scheme. Here we have the representations 

i^ MS >/ + ^ Q 2 , m 2 ) = FS ACT (n f , A, Q 2 , m 2 ) 

-C MP (»/. ^ Q 2 ' ™ 2 ) + *iT("/ + 1. ^ Q 2 ' ™ 2 ) ■ (n.2) 

The pieces in this formulae represent first the results in FOPT, given by 

Q 2 



Q 



lit 



+a 2 (n f , , 2 ) IE S) ® C' (2) ^) 

U=i ^ 



+ft L °(n f ^ 2 )®H?f\9l% 



+al(n f , , 2 ) E 4f£*(nf, S) ® C^' 2 ' (*> ^) ■ (H-3) 



The next pieces are the structure functions FfQ YMP (nf) which can be obtained from 
Fj E g ACT (n/) by replacing all exact heavy quark coefficient functions H i k and L i k (k = q,g) 
by their asymptotic analogues which are defined by 

H^ MP = Jim 2 ^ fc , L™= Jim L hk . (II.4) 

Q 2 >m 2 Q 2 >m. 2 

Finally the structure functions F i P g F (n/ + 1) which are very often called the ZM-VFNS 
representations are defined by 



F™ F (n f + l,A,Q*,m*)=ei 



f^h>f + m 2 )C (0) 



+/r LO (-/+i^ 2 )«c ) (f)} 

+« 2 (-/ + 1, ^ {/&,(«/ + 1, /i 2 ) ® (C (2) (^/ + hp 

^? (a, (^))+EOn/ + l,/' a ))®C W (§) 
r J i=i p 
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+# LO (n / + l,ji 

r, 

,2/ I -i 2\ 2fLO/ 2\ x-x 7- HARD,ASYMP,NS,(2) / A V" 

+a s (n/ + l,/i )}^e k f k+ - k {n f ,n ) ® L i( , V A ' ^2" 



fc=i 



(II.5) 



In all these results the heavy quark Q refers to the bottom quark and the other parameters are 
defined above. For simplicity we will refer to these structure functions as Fp^ ACT (rif = 4), 
Fff YMF (nf = 4), and Ff^ F [rif = 5) respectively, which indicates that the first two structure 
functions depend on four-flavor densities and the last one depends on five-flavor densities. 
The parton densities f k in the above formulae are represented in leading order (LO), next- 
to- leading order (NLO) and next-to- next-to-leading order (NNLO). The NNLO case refers 
to the boundary conditions imposed in |J since the three-loop splitting functions are not 
known yet. These parton densities have been constructed in H starting from a three-flavor 
parametrization in 0. The multiplication of the densities with the heavy and light parton 
coefficient functions is done in such a way that the perturbation series is strictly truncated 
at order a 2 . This is necessary to avoid scheme dependent terms which would otherwise arise 
beyond order a 2 . Therefore the following requirement is satisfied 

^ SN (% = 5)=^ MSN (% = 5 )=^ XAC >/ = 4 ) fOT Q 2 <™ 2 - (H.6) 

Since /<g(m 2 ) NNLO 7^ (see @) this condition can be only satisfied when we truncate the 
perturbation series at the same order. Furthermore because of Eq. ( |ll.4| ) and the property 

^CSN,(l)/ . 1 Q 2 n {l)f . Q 2 N /tt rj\ 

Jim C i>k u (n f + 1, — , — ) = CH(n f + 1, — ) , (II.7) 

we have the asymptotic relation 

lim F™ m (n f + l,A,Q 2 ,m 2 ) = lim F™*(n f + 1, A, Q 2 , m 2 ) 

Q 2 >m 2 ' Q 2 >m 2 

= F™ F (n f + l,A,Q 2 ,m 2 ). (11.8) 

At first sight the form of the expression for -F i B Q ISN in Eq. ( |II.2| ) looks quite different from the 
one presented for -F/q N in Eq. flll.lp . However this is not true. Using the mass factorization 
relations for the asymptotic heavy quark coefficient functions in |2| one can cast -Fj B g ISN into 
the same form as presented for Fpg^ where all quark coefficient functions of the type C 4 C q N 
are replaced by their light quark analogues C^ q appearing in Eq. ( |11.5| ). This replacement 
also applies to the C^q N occurring in the mass factorization relations for C^ 8 N,s and C^j SN ' PS 
presented in ||. Therefore the difference between the CSN and BMSN schemes can be 
attributed to the powers (m 2 /Q 2 y showing up in C^q N but absent in C iA . This effect is only 
noticeable in the threshold region where Q 2 ~ m 2 as we will show below. 

The heavy quark coefficient functions Cff N , H i k) L ijk (k = Q, q, g) and the light partonic 
coefficient functions C^k (k = q, g) can be traced back to the following processes 
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r CSN,S,(l) „S,(1) . * 



^CSN,S,(2) rrS,(2) * . 
Ci,g > H i,g -7+9 



^CSN,PS,(2) „PS,(2) ...*,„ 



Q + Q [6] (CSN), § (EXACT), 
ill (ASYMP) 



Q + Q + 2 1 (CSN), § (EXACT), 
|TJ (ASYMP) 



Q + Q + 9 Bethe-Heitler reaction 



-HARD,NS,(2) 



§ (CSN), (EXACT), Q (ASYMP) 

■y* + q-^Q + Q + q Compton reaction 

f§ (EXACT and ASYMP) 



^CSN.NS.CO) rrNS,(0) . 



: 7 * + Q -> Q 



c. 



CSN,NS,(1) rrNS,(l) . * 



i,<5 



' U i,Q 

,NS,(0) 
<? 

NS,(1) 



:7 * + g_Q + ^ [12 



Z,<J 



^NS,(2) * . 

Ci, q 7 +Q 



^NS,(2) * . 
Ci,q 7 +Q 



/>PS,(2) * . 
Kg : 7 +Q 



CS,(1) * . 
i,s -7 +9 



7 + g -> g 



: 7 * + g^g + S HI 



g + 2 + 3 Pj 



g + g + g 



g + g + g II 



g + g [0 



q + q + g 
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(II.9) 



Behind the reactions we have quoted the references in which the corresponding coefficient 
functions can be found. Note that the heavy quark coefficient functions are mass singular 
when m — > 0. This can be observed immediately when one looks at H^ YMP which behaves 
like ln m (/x 2 /m 2 ) ln n (Q 2 /m 2 ) (see After the logarithms are removed one obtains the 

quantities C^ s which, even though they depend on m, are finite in the limit m — > 0. The 
coefficient function Lf^ RD is finite by itself because as we mentioned above we have imposed 
a lower cut off A = 100 (GeV/c) 2 on the invariant mass of the Q (Q-pair. Finally notice that 
all parton densities, coefficient functions and the running coupling constant are presented 
in the MS-scheme. 

Now we present results for the various structure functions. We are interested in the 
bottom quark structure functions Fp^ N (nf = 5) and i^ B 6 (n^ = 5) for i = 2, L in NNLO 
for the CSN B and BMSN H schemes respectively. In Fig. 1 we have plotted the structure 
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functions F 2 c 6 SN (n/ = 5), F 2 BMSN (n/ = 5), F££ F (n f = 5) and Fg CACT (n/ = 4) in the region 
20 < Q 2 < 10 3 in units of (GeV/c) 2 for x = 0.05. This figure reveals that there is hardly 
any difference between the BMSN and CSN prescriptions. The curves in both prescriptions 
are essentially identical to that for F 2 E ^ ACT (rif = 4). In this region F F ® F (nf = 5) is larger 
than the other results which is expected from the discussion of the bottom quark density 
given in ||. There is still an appreciable difference at the highest plotted Q 2 demonstrating 
that mass effects are important up to very large scales. Notice that for Q 2 < 35 (GeV/c) 2 
F 2 p b DF (nf = 5) becomes negative which means that bottom quark electroproduction cannot 
be described by this quantity anymore. In Fig. 2 we present the same plots for x = 0.005. 
Again one cannot distinguish between F 2 B ^ SN (nf = 5) and F^ (nf = 5) but now both 
are smaller than F 2 Fj b XACT (nf = 4) over the whole Q 2 range. The latter is smaller than 
P^2,b ( n f = 5) in particular for Q 2 > 35 (GeV/c) 2 . Further we want to emphasize that due 
to our careful treatment of the threshold region there is an excellent cancellation between 
F F b DF (n f = 5) and F 2 ASYMP (n/ = 4) so that both F^ N (n f = 5) and F 2 B 6 MSN (n/ = 5) tend to 
F 2 ^ ACT (n / = 4) at Q 2 = m 2 . At large Q 2 we have a cancellation between F 2 A ^ YMP (n f = 4) 
and F^ ACT (n f = 4) so that both F 2 c 6 SN (n/ = 5) and F^ SN (n f = 5) slowly tend to 
F 2 prjF (nf = 5). They are only identical at extremely large Q 2 demonstrating that mass 
effects are still important over a wide range in x and Q 2 . 

In Fig. 3 we show similar plots as in Fig.l for the bottom quark longitudinal structure 
functions. Here we observe a small difference between the plots for F^ N (nj = 5) and 
F^ SN (n f = 5) in the region 20 < Q 2 < 10 3 (GeV/c) 2 . Furthermore F[$ F (n f = 5) is 
considerably larger than the other three structure functions, which differs from the behavior 
seen in Fig.l. This can be mainly attributed to the gluon density which plays a more 
prominant role in than in F 2 ^. For x = 0.005 (see Fig.4) the small difference between 
the BMSN and the CSN descriptions becomes more conspicuous for low Q 2 . 

In Figs.5 and 6 we make a comparison between the NLO and the NNLO structure 
functions F^ N (n f = 5) and F^ SN (n f = 5). Both the CSN and and BMSN descriptions 
lead to the same results in both NLO and NNLO. However while going from NLO to NNLO 
the the structure functions F 2 § N (rif = 5) and F^^ (nf = 5) increase a little bit. The 
differences in the case of x = 0.005 in Fig. 6 are even smaller than those observed for x = 0.05 
in Fig. 5. The same comparison between NLO and NNLO results is made for the longitudinal 
structure functions in Figs. 7 and 8. Here the differences between NLO and NNLO cases 
are much larger than in the case of F 2 ^ in Figs. 5, 6. In NLO both F^^ (nf = 5) and 
F^f N (nf = 5) are smaller than the NNLO results. 

Previous results for F E b XACT (a;, Q 2 , m 2 ) and have been presented in Figs. 20a,20b in |J 
for a now obsolete set of parton densities, so the values quoted there are too small. To 
show these changes we add in Figs. 9 and 10 plots for the x dependence of F 2 c fe SN (n/ = 5), 
F 2 B 6 MSN (n/ = 5), F 2 p fc DF (n / = 5) and F 2 E 6 XACT (n / = 4) at fixed Q 2 = 30 and Q 2 = 100 in units 
of (GeV/c) 2 respectively. Finally we also show in Figs. 11 and 12 plots for the x dependence 
of ifJN( n/ = 5 ) ; F™ SN (n f = 5), Fl$ F (n f = 5) and F F f ACT (n f = 4) at fixed Q 2 = 30 and 
Q 2 = 100 in units of (GeV/c) 2 respectively. Note that there are also some recent results for 
F 2 ,b in H in the TR scheme and in jL4j] for FOPT. 

The plots for F^^ N (nj = 5) in Figs. 7,8 do not show a negative region at small Q 2 which 
could have been expected by analogy with the results for F^^(nf = 4) in |J. In the case 
of bottom quark production the negative regions do occur but at even smaller values of x. 
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In Figs. 13, 14 we show the same plots as in Figs. 1,3 respectively but for x = 5 x 10~ 5 . 
Now the structure function F^f N is negative in the region Q 2 « 30 (GeV/c) 2 . In Figs. 
15,16 we show the same plots as in Figs. 5,7 respectively but for x = 5 x 10~ 5 . Fig. 16 
shows that the longitudinal structure functions for the case of bottom production also have 
negative regions at small Q 2 values in both NLO and NNLO. This phenomenon also occurs 
for the charm quark structure functions in ||. In the NLO case this arises because the term 
/S.NLO( n/ + 1? ^ g, c csn,s,(i) (g2/ m 2 ; q2/ // 2\ in Eq is negative due to the definition 



of the gluon coefficient function in the CSN scheme (see Eq.(2.19)) in ||) which is given by 



^csn,s,(i)/Q 2 Q 2 
W> 2 



H 



ni 



,S,(1)//^ n r CSN,NS,(0)/y \ 



with C 



CSN,NS,(0) 
L,Q 



Am 



(11.10) 



where Aq'^ denotes the one- loop operator matrix element computed in [[□]]. Notice that 

s (1) 

the latter and the lowest order exact coefficient function H Lg are always positive. Because 
of the minus sign in Eq. ( |1L10| ) it appears that the coefficient function £CSN,s,(i) become 
negative in particular at low Q 2 values. In the NNLO case one obtains more negative 
contributions due to the term fn^k°(nf + 1,// 2 )C£^' NS in formula ( [11. 1| ). It turns out 
that fQ+Q°(nf + l,x, /i 2 ) is negative at small x and /i 2 = Q 2 > m 2 . Notice that at the 
latter scale /q+q(?V + /i 2 ) and /q+^Ov + l,x,/i 2 ) are very small because they vanish 
at /i = m in contrast to /q+^°(^/ + l,a;, /i 2 ). The behavior of the structure function above 
is characteristic of the CSN scheme since it does not appear in the case of BMSN. This is 
because in the latter scheme the longitudinal coefficient function, represented by £CSn,ns,(o)^ 
is identical to zero so that the zeroth order contribution to (rif = 5) vanishes and the 

first order correction is given by £BMSN,s,(i) = H^ g . The latter leads to a positive structure 
function over the whole kinematical region. To further demonstrate this point we plot in 
Fig. 17 pieces of the NLO result 



FE™(n f + l,Q\m 2 ) = el 



/S?(n/ + l,A* a )CS 



CSN,NS 



(o)/Q 



=5) 



+0,(71, + 1, /1 2 ) \ fl${n f + 1, /x 2 ) <g> C 



7 CSN,NS 
i,6 



(1) / Q 2 Q 2 



+ /r( % + i,^ )s c5 NA(1, (S-f)}' 



(11.11) 



The sum of the b-quark contributions, labelled Terml, is always positive. The gluonic con- 
tribution, labelled Term2, is clearly negative over a wide range in Q 2 and large enough that 
the sum F^f N is also negative for 30 < Q 2 < 150 GeV 2 , as in Fig.16. This behavior is due to 
the gluonic coefficient function £CSN,s,(i) w hich is explained below Eq. (|H.10|) . However the 
magnitude of the gluonic contribution depends on the choice of the gluon density. If we use 
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an NLO gluon density in the order a s contribution to the structure function in Eq. (|H.11|) 
rather than a LO gluon density then the sum of the first two terms is unchanged but the 
gluonic part is now smaller in magnitude. These contributions are shown in Fig. 18 where 
now |Term2| < Terml so that the total result for -Ff^ N is everywhere positive. However this 
procedure violates our prescription for the computation of the structure functions in both 
the CSN and the BMSN schemes. In this prescription the LO densities are multiplied by the 
highest order coefficient function whereas the NLO densities are combined with lower order 
coefficient functions (see formulae ( |IL1|) , (|II.3|) and (|II.5|) ). In this way the perturbation 
series is truncated up to the order we want to compute the structure functions. Hence we 
avoid terms, arising beyond that order, which introduce a scheme dependence and spoils the 
threshold behavior (see f|). The latter happens if one follows the usual procedure where 
one multiplies the highest order densities by the highest order coefficient functions. The 
difference between the usual procedure and our prescription is not only shown by our parton 
density set but is also observed for other sets presented in the literature. Examples are 
recent sets such as MRST98 (with m b = 4.3 GeV, m c = 1.35 GeV), MRST99 p| (with 
m b = 4.3 GeV, m c = 1.43 GeV), and CTEQ5 (with m b = 4.5 GeV, m c = 1.3 GeV). 
Note that the MRST99 set does not provide LO densities. Using their NLO densities they 
yield positive values for the Q 2 dependence of F^ b at x = 5 x 10~ 5 . There are both LO 
(CTEQ5L) and NLO (CTEQ5M) densities in the CTEQ5 set and we have checked that, for 
the same x, Q 2 values F^f N is positive with purely NLO densities but has a negative region 
when the LO and NLO densities are used according to our prescription. The observations 
made above leads to the conclusion that the 4m 2 /Q 2 term in the non-singlet CSN longitudi- 
nal coefficient function in Eq. ( [11. 10| ) leads to a negative gluonic coefficient function. When 
the latter is used together with the latest LO and NLO parton density sets F^f^(nf = 5) 
becomes negative in the low Q 2 -region at small x. We speculate that the CSN scheme would 
always yield positive structure functions if we could use parton density sets which fitted 
data with convolutions of LO densities with 0(a s ) coefficient functions and NLO densities 
with zeroth order coefficient functions. Unfortunately such densities are not available. 

To summarize the main points we have implemented two variable flavor number schemes 
for bottom quark electroproduction in NNLO and compared them with NLO FOPT results. 
The schemes differ in the way mass factorization is implemented. In the CSN scheme this 
is done with respect to the full heavy and light quark structure functions at finite Q 2 . In 
the BMSN scheme the mass factorization is only applied to the coefficient functions in the 
large Q 2 limit. Both schemes require four-flavor and five-flavor parton densities which satisfy 
discontinuous NNLO matching conditions at a scale /i = m b . We have constructed these 
densities using our own evolution code M. The schemes also require matching conditions 
on the coefficient functions which are implemented in this paper. Note that we have also 
removed the dangerous terms in \v?{Q 2 /m 2 ) from the Compton contributions so that both 
Fp b SN (nf = 5) and F^^ (rif = 5) are collinear safe. As in we have done this in a way 
which is consistent with our study of inclusive quantities by implementing a cut A on the 
mass of the b — b pair. We stress that any ZM-VFNS bottom quark density description of 
Fif, must use collinear safe definitions. This is not required in the fixed order perturbation 
theory approach given by Ff b XACT (rif = 4) in JI| for moderate Q 2 -values. 

Finally we made a careful analysis of the threshold behaviors of Fj C b SN (n/ = 5) and 
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Fj B fe MSN (n/ = 5). In order to achieve the required cancellations at the scale \i = nib so that 
they both become equal to Ff^ {rif = 4) one must be very careful to combine terms 
with the same order in the expansion in a s . The approximation we made in this paper, of 
using NLO splitting functions in place of NNLO splitting functions, was sufficient for our 
purposes. We successfully implemented the required cancellations near threshold and the 
corresponding limits at large scales came out naturally. Inconsistent sets of parton densities 
automatically spoil these cancelations. Since there are only minor differences between the 
CSN, BMSN and NLO FOPT predictions it is clear that the use of variable flavor number 
schemes for bottom quark production is not required for the analysis of HERA data. However 
the ZM-VFNS description is clearly inadequate at small Q 2 . 
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Figure Captions 

Fig. 1. The bottom quark structure functions ^2& CACT ( 7 V = 4) (solid line) F^^^rif = 5), 
(dot-dashed line) F^ SN (n f = 5), (dashed line) and F$ F (n f = 5), (dotted line) in 
NNLO for x = 0.05 plotted as functions of Q 2 . 

Fig. 2. Same as in Fig. 1 but now for x = 0.005. 

Fig. 3. The bottom quark structure functions F F f ACT (nf = 4) (solid line) F^f N (nf = 5), 
(dot-dashed line) F^ SN (n f = 5), (dashed line) and Fl^ F (n f = 5), (dotted line) in 
NNLO for x = 0.05 plotted as functions of Q 2 . 

Fig. 4. Same as in Fig. 3 but now for x = 0.005. 

Fig. 5. The bottom quark structure functions F^ SN (n f = 5) in NLO (solid line), NNLO 
(dotted line) and F^ N (n f = 5) in NLO (dashed line), NNLO (dot-dashed line) for 
x = 0.05 plotted as functions of Q 2 . 

Fig. 6. Same as in Fig. 5 but now for x = 0.005. 

Fig. 7. The bottom quark structure functions F^ SN (nf = 5) in NLO (solid line), NNLO 
(dotted line) and Ff| N (n/ = 5) in NLO (dashed line), NNLO (dot-dashed line) for 
x = 0.05 plotted as functions of Q 2 . 

Fig. 8. Same as in Fig. 7 but now for x = 0.005. 

Fig. 9. The bottom quark structure functions F^ ACT (nj = 4) (solid line) F^^irif = 5), 
(dot-dashed line) F^ SN (n f = 5), (dashed line) and F^ F {n f = 5), (dotted line) in 
NNLO for Q 2 = 30 (GeV/c) 2 plotted as functions of x. ' 

Fig. 10. Same as in Fig. 9 but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 11. The bottom quark structure functions F^f ACT (nf = 4) (solid line) F^(nf = 5), 
(dot-dashed line) F^ SN (n f = 5), (dashed line) and F^ijif = 5), (dotted line) in 
NNLO for Q 2 = 30 (GeV/c) 2 plotted as functions of x. 

Fig. 12. Same as in Fig. 11 but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 13. Same as in Fig. 1 but now for x = 5 x 10~ 5 . 

Fig. 14. Same as in Fig. 3 but now for x = 5 x 10~ 5 . 

Fig. 15. Same as in Fig. 5 but now for x = 5 x 10~ 5 . 

Fig. 16. Same as in Fig. 7 but now for x = 5 x 10~ 5 . 

Fig. 17. The bottom quark structure function F^ ACT (nf = 4) (solid line) F^^ N (n/ = 5), 
(dot-dashed line) together with the NLO charm density piece Terml (dotted line) and 
the LO gluon density piece Term2 (dashed line), see text, for x = 5 x 10^ 5 plotted as 
functions of Q 2 . 
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Fig. 18. The bottom quark structure function Ff^ ACT {rif = 4) (solid line) F^f N (nf = 5), 
(dot-dashed line) together with the NLO charm density piece Terml (dotted line) and 
the NLO gluon density piece Term2 (dashed line), see text, for x = 5 x 10~ 5 plotted 
as functions of Q 2 . 
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